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Universality in the synchronization of weighted random networks 
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Realistic networks display not only a complex topological structure, but also a heterogeneous 
distribution of weights in the connection strengths. Here we study synchronization in weighted 
complex networks and show that the synchronizability of random networks with large minimum 
degree is determined by two leading parameters: the mean degree and the heterogeneity of the 
distribution of node's intensity, where the intensity of a node, defined as the total strength of input 
connections, is a natural combination of topology and weights. Our results provide a possibility for 
the control of synchronization in complex networks by the manipulation of few parameters. 
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Over the past few years, the analysis of complex sys- 
tems from the viewpoint of networks has become an im- 
portant interdisciplinary issue jj. It has been shown 
that physical and dynamical processes, such as cascad- 
ing failures 2|, epidemic spreadings 3], and network syn- 
chronization _4, 5 , ii, Xi iSj i2(| J are strongly influenced by 
the structure of the underlying network. Previous work 
on synchronization has focused mainly on the influence 
of the topology of the connections by assuming that the 
coupling strength is uniform. However, synchronization 
is influenced not only by the topology, but also by the 
strength of the connections ^9j . Most complex networks 
where synchronization is relevant are indeed weighted. 
Examples include brain networks 10], networks of cou- 
pled populations in the synchronization of epidemic out- 
breaks [ll|, and technological networks whose function- 
ing relies on the synchronization of interacting units |l2l | . 
The distribution of connection weights in real networks 
is often highly heterogeneous [l^- The study of syn- 
chronization in weighted networks is thus of substantial 
interest. 

In this Letter, we address this question in random net- 
works with weighted coupling schemes motivated by real 
networks. Our main result is the uncovering of a univer- 
sal formula that describes with good approximation the 
synchronizability of identical oscillators solely in terms of 
the mean degree and the heterogeneity of the node's in- 
tensity, irrespective of the degree distribution and other 
topological properties. The intensity of a node, defined 
as the sum of the strengths of all input connections of 
that node, incorporates both topological and weighted 
properties and raises as a very important parameter con- 
trolling the synchronizability. In particular, it follows 
that the synchronizability is significantly enhanced when 
the heterogeneity of the node's intensities is reduced. 

The dynamics of a general weighted network of N cou- 



pled identical oscillators is described by 



^N 



X, = F(x,) + ^E LiW/y^.,[H(x,) - H(xO], (1) 



^N 



F(xO-aELiG,,H(x,), i = l,...,N, (2) 



where F = F(x) governs the dynamics of each individ- 
ual oscillator, H = H(x) is the output function, and a 
is the overall coupling strength. Here G = [Gij) is the 
coupling matrix combining both topology [adjacency ma- 
trix A = [Aij)] and weights [weight matrix W = {Wij), 

Wij > 0]: G^j = SijSi-WijAij, where St = Y.'j'^i ^ij^ij 
denotes the intensity of node i. The rows of G have zero 
sum, and this ensures that the completely synchronized 
state {xi = s,V i |s = F(s)} is an invariant manifold of 
Eq. |(2Jl. In this work, we focus on the class of weighted 
networks where G is diagonalizable and has real eigenval- 
ues. As it will be shown shortly, in this case the synchro- 
nizability of the networks can be characterized by the 
properties of the eigenvalues, without referring to spe- 
cific forms of F and H. This applies in particular to the 
important class of networks where G can be written as 
G = BC for B a nonsingular diagonal matrix and C a 
symmetric, zero row-sum matrix. We assume that ma- 
trix A is binary and symmetric, and that the connection 
weights and asymmetries are incorporated into W. 

The linear stability of the synchronized states can 
be studied by diagonalizing the variational equations of 
Eq. lO into N blocks of the form ^ 



ii = [DFis) ^ a\iDIlis)]Ci, 1 = 1, 



.N, 



(3) 



which are different only by A;, 
ordered as = Ai < A2 • • • < 
G are nonnegative because Gu 



the Ith eigenvalue of G, 
Xn- The eigenvalues of 
= - J2j^^ Gij > and 



Ai = because J2i ^v — for all i. The stability of s is 
determined by A; and the master stability function Il4|, 
i.e., the largest Lyapunov exponent A of the generic vari- 
ational equation ^ = [£'F(s) — eDH(s)]^. For many os- 
cillatory dynamical systems |y, |lj|, A is negative in a 



single, finite interval ei < e < £2, where the thresholds 
ei and €2 are determined only by F, H, and s. The net- 
work is thus synchronizable for some a iff the condition 
ei < crA; < 62 is satisfied so that A((tA;) < for all I > 2. 
This is equivalent to the condition 



R = Xn/>^2 < £2/61, 



(4) 



where the eigenratio R depends only on the network 
structure (G), and £2/61 depends only on the dynamics. 
From these, it follows that the smaller the eigenratio R 
the more synchronizable the network and vice versa [6j. 
Another measure of synchronizability is the cost C in- 
volved in the couplings of the network |9| . When Eq. Q 
is satisfied, the synchronized state is linearly stable for 
f > finin = £1/^2- The cost C is the total input strength 
of the connections of all nodes at the synchronization 
threshold: C = cr,nin Z^ij ^^u^y" = (■ 
normalized cost 



^N 



EtiS^- The 



Co = C/iNei) = n/X2, 



(5) 



^N 



where fl = X]i=i ^if-^i does not depend on the dynam- 
ics (F, H, and s) and can be used as a complementary 
parameter of synchronizability 9]. Here we character- 
ize the synchronizability of the networks using both the 
eigenratio R and the cost Cq. 

Previous work has obtained bounds for the eigenvalues 
of unweighted networks {Wij = 1) (see [l5| for a review). 
Such bounds, however, are not tight and may provide 
limited information about the actual synchronizability 
of complex networks. Here we aim at obtaining a more 
quantitative approximation of the synchronizability for 
a class of weighted random networks with real spectra, 
which includes unweighted networks as a special case. 
Our analysis is based on the combination of a mean field 
approximation and new graph spectral results. 

First, in random networks with fcmin 3> 1, close to a 
synchronized state Eq. |^ can be approximated as 

i, = F(x,) + (7(5,/fc,0Ef=i^y[H(x,) - H(x,)]. (6) 

The reason is that each oscillator j receives sig- 
nals from a large and random sample of other os- 
cillators in the network and Xj is not affected di- 
rectly by the individual output weights Wij. Conse- 
quently, we may assume that Wij and H(xj) are sta- 
tistically uncorrelated and that J2i=i^ij^ij'^i^j) ~ 



•.N 



•.N 



(i/fc»)E;=iW^u-4^jE;=i^.jH(x,) = 5,H, if fc, » 1 

[l6|. Here H,; — {1/ki) Ei=i ^iiH(xj) is the local mean 
field. 

Now, if the network is sufficiently random, the local 
mean field H^ can be approximated by the global mean 
field of the network, H^ « H = (1/7V) ^^^^ H(xj). 
Moreover, close to the synchronized state s, we may as- 
sume Hi w H(s), and the system is approximated as 



indicating that the oscillators are decoupled and forced 
by a common oscillator s = F(s) with forcing strength 
proportional to the intensity Si. The variational equa- 
tions of Eq. {7|) have the same form of Eq. Q, except 
that A; is replaced by Si. If there exists some a satis- 
fying ei < aSi < 62 for all I, then all the oscillators are 
synchronizable by the common driving H(s), correspond- 
ing to a complete synchronization of the whole network. 
These observations suggest that the eigenratio and the 
cost can be approximated as 



H ~ >-'max/»-'minj ^0 ~ ' '/'-'n 



(8) 



maximum 



where ^min, 'S'max, and Q are the minimum, 
and mean intensities, respectively. 

Next we present tight bounds for the above approx- 
imation. Eq. lO means that the coupling matrix G is 
replaced by the new matrix G" = (G°), with G° — 
f-iSijki - Aij). G" can be written as G° ^ SG = 
SD-\D - A), where S = (SijS^) and D = {6ijh) are 
the diagonal matrices of intensities and degrees, respec- 
tively, and G is the normalized Laplacian matrix |l7| . 
Importantly, now the contributions from the topology 
and weight structure are separated and accounted by G 
and S, respectively. We can show that the largest and 
smallest nonzero eigenvalues of matrix G" are bounded 
by the eigenvalues /i/ of G as 

SminfJ-2C < A2 < OminC , >5max !i ^N < Omax/^Af; (9) 

where c and c' can be approximated by 1 for most large 
complex networks of interest, such as realistic scale- free 
networks (SFNs). The proof is involved and long, and 
the details will be presented elsewhere. If the network 
is sufficiently random, the spectrum of G tends to the 
semicircle law for large networks with arbitrary expected 
degrees !l3|, provided that fc,„i„ 3> \/l<, and max{l — 
^^2,^^-1} = [l + o(l)]^ for fcmin > Vi^ln^A^, where 
K is the mean degree. From these, it follows that 



M2 



1 - 2/Vk, fiN^l + 2/Vk, 



(10) 



which we find to provide a good approximation under the 
weaker condition fcmin ^ 1 , regardless of the degree dis- 
tribution. From Eqs. © and (|10|l . we have the following 
approximations for the bounds of R and Cq: 

Jmax ^ r, ^ '^max 1 + ^/ V-iV /i i \ 

S -fi < — -^=, (ttj 



»^inii] 



< Go < 



S'min 1~2/Vk' 

n 1 

S'min 1-2/Vk' 



(12) 



X, =F(x,) + fT5,[H(s)-H(x,)], 



1, 



.N, (7) 



For the case of unweighted networks {Si = ki, ^ = K), 
the bounds in Eq. Hll|l are much tighter than those re- 
viewed in Ref. \l^ . 

The bounds in Eqs. Hll|) and (|12|) show that the con- 
tribution of the network topology is mainly accounted 



by the mean degree K . Therefore, for a given K, the 
synchronizabiUty of random networks with large fcmin is 
expected to be well approximated by the following uni- 
versal formula: 



R = An^, Co = Ac-^ 



On 



On 



(13) 



where the pre-factors Aji and Ac are expected to be 
close to 1. In the case of matrix G" with uniform inten- 
sity [Si = 1 Vz), they are given by the upper bounds, 
^« = T^ -d Ac = ^^^, and A^ -. 1 and 
Ac — > 1 in the limit K ^ oo. Formula Ijl^f) is con- 
sistent with the approximation in Eq. and indicates 
that the synchronizabiUty of these networks is primarily 
determined by the heterogeneity of the intensities. 

Our numerical results on various weighted and un- 
weighted networks have confirmed this universal formula. 
First we consider the following weighted coupling scheme: 



W,, 



Di/ ki 



(14) 



in which the intensities Si follow an arbitrary distribu- 
tion not necessarily correlated with the degrees. In this 
case, Eq. ((HJ and Eq. QJ are identical and C^ — G. 
This weighted coupling scheme includes many previously 
studied systems as special cases. If Si = ki Vi, it cor- 
responds to the widely studied case of unweighed net- 
works la, IE lal- I'^ ttis case of fully uniform intensity 
{Si = 1 Vi), it accommodates a number of previous stud- 
ies about synchronization of coupled maps |j, |3|. The 
weighted scheme studied in |2| , Wij = kf , is another spe- 
cial case of Eq. (dj) where Si = k]^^ . 

We have applied the weighted scheme to various net- 
work models: (i) Growing SFNs with aging [iSj. Start- 
ing with 2m + 1 fully connected nodes, at each time step 
we connect a new node to m existing nodes according 
to the probability 11^ ~ hT"" , where n is the age of 
the node. The minimum degree is then fcmin = fn and 
the mean degree is K — 2m. For the aging exponent 
— oo < a < 0, this growing rule generates SFNs with 
a power-law tail P{k) ~ fc~^ and scaling exponent in 
the interval 2 < 7 < 3 [ll|, as in most real SFNs. For 
a = 0, we recover the usual BA model |]j|, which has 
^ — 3. (a) Random SFNs |20j. Each node is assigned to 
have a number ki > fcmin of "half-links" according to the 
distribution P{k) ~ k^'^ . The network is generated by 
randomly connecting these half-links to form links, pro- 
hibiting self- and repeated links. (Hi) K-regular random 
networks. Each node is randomly connected to K other 
nodes. 

We present results for two different distributions of in- 
tensity Si which are uncorrelated with the distribution of 
degree kf. (1) a uniform distribution in [^mim "Smax]; and 
(2) a power-law distribution, P{S) ~ S~^ , S > 5'min, 
where S'min is a positive number. Consistently with the 




FIG. 1: (a) _R as a function of S'max/S'niin and (b) Co as a func- 
tion of il/Smin, averaged over 50 realizations of the networks. 
Filled symbols: uniform distribution of Si G [5'min, Smax]. 
Open symbols: power-law distribution of Si, P{S) ~ S" 
for 2.5 < r < 10. Different symbols are for networks with dif- 
ferent topologies: BA growing SFNs (o), growing SFNs with 
aging exponent a — —3 (□), random SFNs with 7 = 3 (o), 
and A-reguIar random networks (A). The number of nodes 
is A = 2'^'' and the mean degree is A = 20. Insets of (a) and 
(b): Ar and Ac as functions of A for Smax/Smin = 1 (o), 2 
(D), 10 (A), and 100 (*), obtained with uniform distribution 
of Si in A-regular networks. The dashed lines are the bounds. 
Solid lines in (a) and (b): Eqs. ifT!^ with Ar = Ac = 1. 



prediction of the universal formula, if fcmin 3> 1, the 
eigenratio R collapses into a single curve for a given K 
when plotted as a function of S'max/5'min [Fig. 1(a)], irre- 
spective of the distributions of ki and Si . The same hap- 
pens for the cost Co as a function of O/S'min [Fig- 1(b)]. 
The behavior of the fitting parameters Aj^ and Ac is 
shown in the insets of Fig. 1. For uniform intensity, 
they are very close to the upper bounds. They approach 
very quickly 1 when the intensities become more het- 
erogeneous (S'max/S'min > 3). Therefore, Eqs. ((T^ with 
Aji = Ac — 1 [Fig. 1, solid lines] provide a good approx- 
imation of the synchronizability for any large K if the 
intensities are not very homogeneous. 

In more realistic networks, including scientific collabo- 
ration networ ks 1 131] , metabolic networks |21| , and airport 
networks il3|, |2l| , it has been shown that the weight Wij 
of a connection between nodes i and j is strongly corre- 
lated with the product of the corresponding degrees as 
(Wij) ^ (kikj)^ . Here 6 depends on the specific net- 
work under study. Motivated by these observations, we 
analyze the weighted coupling 



W.,, = ihk.f 



(15) 



where the weights are defined for the connections of a 
given network topology and is a tunable parameter. 9 
controls the heterogeneity of the intensity Si and the cor- 
relation between Si and ki, since Si — k]^ {k^j)i, where 
{k^j)i = (l/ki) J2 ^f is approximately constant for fc^ ^ 1 
when the degree correlations can be neglected. Variations 
of 9 have significant impact on the synchronizability of 
SFNs [Fig. 2]. However, as shown in the insets of Fig. 2 
for various networks and 9 values, R and Cq collapse 
again to the universal curves when regarded as functions 
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FIG. 2: (a) Eigenratio R and (b) cost Co as functions of 6 for 
the BA growing SFNs (o) , growing SFNs with aging exponent 
a = —3 (□) and random SFNs with 7 = 3 (A). Each symbol 
is an average over 50 reahzations of the networks with K — 20 
and A*' — 2^". Inset of (a): the same data for i? as a function 
of Smax/Smin- lusct of (b): the Same data for Co as a function 
of n/Smin- Sohd lines: Eqs. (E3J with Ar = Ac ^ 1. 
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FIG. 3: (a) i? as a function of the system size N for un- 
weighted networks {9 = 0). (b) The same data for i? as a 
function of Smax/Smin. Solid line: Eq. CHJ with Ar = 1. The 
other parameters and legends are the same as in Fig. 2(a) 



of 5'niax/S'min and fl/Snjin, respectively. 

We emphasize that these results also hold for un- 
weighted random networks. As shown in Fig. 3 for un- 
weighted SFNs, the eigenratio _R as a function of the 
network size N collapses to a single universal curve when 
plotted against S'max/S'min •^■ 

In addition to heterogeneous degrees and weights, 
many real networks also display hig h clustering ' 23| and 
nontrivial correlation of degrees |24| . The clustering and 
degree correlations are negligible in random SFNs and 
-fC-regular random networks, but we find that they are 
significant in growing SFNs for a < 0. The results in 
Figs. 1-3 show that the universality holds without signif- 
icant dependences on these topological properties. 

Eqs. H13|l also provide a meaningful approximation for 
networks which are not fully random. For example, con- 
sider small- world networks |23| where a regular ring of N 
(= 2^") nodes, each connected to K (= 20) nearest neigh- 
bors, is rewired with a probability p for each link. We 
find that R and Co collapse to the universal curves even 
when the networks are dominated by local connections, 
e.g. for p — 0.3, if the intensities are very heterogeneous 
(S'max/S'min > 10). For networks with fcmin ~ 1, the syn- 
chronizability is still strongly dependent on Smax/Smin 
and if, although it shows additional dependences on the 
details of the distributions of Si and hi and on other 
topological properties. 



In summary, we have shown that the synchronizabil- 
ity of sufficiently random networks with minimum de- 
gree fcmin ^ 1 is universally dominantly determined by 
the mean degree K and the heterogeneity of the inten- 
sities Si. This universality applies to a general class of 
large networks where the heterogeneity of Si is due to ei- 
ther the distribution of degrees, as in unweighted SFNs, 
or the distribution of connection weights, as in weighted 
iiT-regular networks, or a combination of both, as ex- 
pected in most realistic networks, such as in the airport 
network |l3l | , which underlies the synchronization of epi- 
demic outbreaks jlll |. In particular, formula (|13|l explains 
why synchronizability is improved when the heterogene- 
ity of Si is reduced, which can be useful for network de- 
sign and control of synchronization. 
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